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ABSTRACT 

This paper is the second part of a work investigating the properties of dusty tori in 
\ active galactic nuclei (AGN) by means of multi-component spectral energy distribution 

(SED) fitting. It focuses on low luminosity, low redshift (z ^ 0.25) AGN selected 
among emission line galaxies in the overlapping regions between SWIRE and SDSS 
Data Release 4 as well as X-ray, radio and mid-infrared selected type 2 AGN samples 
from the literature. The available multi-band photometry covers the spectral range 
' from the w-band up to 160um. Using a standard \ 2 minimisation, the observed SED 

of each object is fit to a set of multi-component models comprising a stellar component, 
^3 a high optical depth (V9.7 > 1.0) torus and cold emission from a starburst (SB). The 

torus components assigned to the majority of the objects were those of the highest 
optical depth of our grid of models (rg.7 = 10.0). The contribution of the various 
components (stars, torus, SB) is reflected in the position of the objects on the IRAC 
ON ! colour diagram, with star-, torus- and starburst-dominated objects occupying specific 

areas of the diagrams and c omposite objects lying in be tween. The comparison of type 
' 1 (as derived from Paper l. lHatziminaoglou et al.ll2008[ ) and type 2 AGN properties is 

broadly consistent with the Unified Scheme. The estimated ratio between type 2 and 
t^. ■ type 1 objects is about 2-2.5:1. The AGN accretion-to-infrared luminosity ratio is an 

' indicator of the obscuration of the AGN since it scales down with the covering factor. 

ON . We find evidence supporting the receding torus paradigm, with the estimated fraction 

of obscured AGN, derived from the distribution of the covering factor, decreasing with 
increasing optical luminosity (AL5100) over four orders of magnitude. The average star 
formation rates are of ~ IOMq/ut for the low-z sample, ~ AOM^/yr for the other type 
' 2 AGN and ~ 115M /yr for the quasars; this result however, might simply reflect 

observational biases, as the quasars under study were one to two orders of magnitude 
more luminous than the various type 2 AGN. For the large majority of objects with 
70 and/or 160 um detections an SB component was needed in order to reproduce the 
data points, implying that the far- infrared emission in AGN arises mostly from star 
formation; moreover, the starburst-to-AGN luminosity ratio shows a slight trend with 
increasing luminosity. 
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1 INTRODUCTION 

Dust is the cornerstone of the active galactic nuclei (AGN) 
Unified Scheme that postulates that the diversity of the ob- 
served properties of AGN are merely a result of the different 
lines of sight with r espect to obscuring material surrounding 
the active nucleus (|Antonuccilll993l : lUrri fc Padovanilll995l ; 
Tadhuntcr 2008). The elements necessary for understanding 
the nature of AGN and the variety of their properties are in 
fact the geometry of the circumnuclear dust and the amount 



of obscuration. The latter is also needed in order to under- 
stand the physical conditions of the dust enshrouded nuclear 
region and to estimate the intrinsic UV-to-optical properties 
of the AGN. Obscuring dusty material should re-radiate in 
the infrared (IR) the fraction of accretion luminosity it ab- 
sorbs, providing thus a direct metric for the study of the 
medium. Indeed all spectroscopically confirmed AGN lying 
in regions covered by the various Spitzer Space Telescope 
(Spitzer) surveys show significant IR emission down to the 
detection limits. 



2 Hatziminaoglou, Fritz & Jarrett 



The conventional wisdom is that dust, distributed in 
a more or less toroidal form around the nucleus, mainly 
consists of silicate and graphite grains, leaving unmistak- 
able signatures in the observed SEDs of AGN. The most 
important dust features are the T~1500 black-body like 
rise of the IR emission at A ^ 1 um, corresponding to 
the sublimation temperature of the graphite grains and the 
~9.7 micron feature in emission/absorption attributed to 
silicates. While silicates in absorption were already long 
ago observed in a variety of type 2 objects, ISO observa- 
tions only provided faint indications of their presence in 
emission in the spectra of type 1, challenging the Unified 
Scheme, as all smooth torus models predicted a feature 
in emission for the unobscured AGN. Clumpy models, on 
the o ther hand, provided a natural attenuation of the fea- 
ture (|Nenkova et all 2002). Observations of tens of objects 
carried out with the IRS spectrograph onboard Spitzer al- 
lowed, how ever, the detection of the silicate feature (first 
reported bv lSiebenmorgen et al.|[2005l and observed by sev- 
eral others since) predicted by smooth tori models in emis- 
sion for type 1 AG N and quasars (e .g. IPier fc Krolik|[l992l : 
iGranato fc Danesdll994l ; iFritz et all 120061 ). Thus far, there 
is no consensus on the dust distribution among astronomers; 
for an overview of the issues involvin g the two appro aches 



. Dullemond fc van Bemmell (|2005l) . lElitzurl ((2008) and 
iNenkova et all (|200Sh . 



Until very recently the number of spectroscopically con- 
firmed type 2 AGN was rather limited; in the last few years 
however, the advent of large-scale multi-wavelength surveys 
has led to a revolution in AGN studies. We are now in a 
position to compile statistically large type 2 AGN samples 
at all redshifts, selected in a multitude of ways. In this work 
we try to characterise the observed spectral energy distri- 
bution (SED) o f type 2 AGN usi ng the smooth torus mod- 
els described in IFritz et all (120061 ) and the methodology de- 
scribed in lHatziminaoglou et alj l|2008l). hereafter Paper 1, 
in which we dealt with the IR properties of quasars. The 
aim of the present work is to model a large sample of type 
2 AGN and, by comparison with the results derived from 
the study of type 1 quasars in Paper 1, test the Unified 
Scheme. The various type 2 AGN samples will be discussed 
in detail in Section[5] Section[3]will briefly describe the var- 
ious model components as well as the basics of the standard 
X 2 minimisation technique and assumptions. The results on 
the individual samples will be presented in Section U while 
the combined results from all samples, including the quasars 
(typical representatives of type 1 AGN) from Paper 1, will 
be shown in Section [S] A final discussion on the future of 
such studies is presented in Section [6] 



2.1 Low-luminosity low-redshift AGN 

The low-luminosity low-redshift AGN sample, hereafter low- 
z sample, is part of an updated version of an emission line 
galaxies sample classified as AGN. The selection was made 
among hundreds of thousands of emission line galaxies from 
the S PSS Data Release 4 (DR4; lAdelman-McCarthv et all 
2006), based on the relative strength of the emission lines 
on a log([Olll]/H/3) - log([Nll]/Ha) diagram, after care- 
ful subt raction of the ste l lar co ntinuum, described in de- 
tail on iKauffmann et al.l (|2003h Q. From the more than 
80,000 proposed AGN, 420 reside in the northern SWIRE 
l|Lonsdale et al1l2003l ; lLonsdale et al 1I2004D fields ELAISN1, 
ELAIS N2 and Lockman. Among these objects, 388 have 
been detected in at least two out of the four IRAC bands, 
with 340 having also MIPS detections at 24 unfl From the 
388 (340 with an additional 24 um detection) objects, 155 
(146) and 88 (81) were additionally detected by MIPS at 70 
and 160 um, respectively. The far-infrared (FIR) measure- 
ments characterize or constrain (with upper limits) the cold 
dust emission. 

Due to the nature of the sample, particular care was 
taken in the construction of the SEDs. We used data 
from the SDSS DR4, 2MASS including the deeper data set 
2MA SS x 6 taken in the Lockman field fecichma n et alj 
2003), as well as IR data from SWIRE. For extended objects 
(based on the SDSS classification) we used Petrosian u, g, r, 
i, z magnitudes, 7" radius magnitudes for 2MASS JHK, and 
aperture 5 fluxes (5.8" radius), the largest publicly available 
IRAC aperture, for the four IRAC bands, to better capture 
the total light arising from the resolved host galaxies. The 
5" 2MASS magnitudes, closer to the 5.8" aperture selected 
for the IRAC bands, were proven to be inadequate when 
all the SEDs were visually inspected as they were clearly 
falling short of both the SDSS and the IRAC fluxes. For 
the point sources we used PSF SDSS magnitudes, default 
(PSF) 2MASS magnitudes and aperture 2 (1.9") IRAC mag- 
nitudes (for a full description of the SWIRE magnitudes 
and their adequat e for point-like and extended sources see 
ISurace et al1l2005l ). The redshift histogram of the 420 ob- 
jects is shown in Fig. [1] with the gray histogram corre- 
sponding to the 388 objects that will be studied in detail. 
Tables [T] and [5] show the SDSS names, coordinates, redshifts, 
emission line properties, and stellar masses (taken from the 
narrow emission line object catalogue), and SWIRE fluxes, 
respectively for the 420 objects. The objects belonging to 
the Subset (S) of 388 entries used for the study are marked 
in the column "S" of Table [T] The objects marked in column 
"AGN" are those that were assigned with an AGN compo- 
nent as a result of the fitting procedure. 



2 THE SAMPLES 

In this work we will be analysing a variety of type 2 AGN 
samples, namely a low redshift, low luminosity AGN sample, 
an X-ray and radio selected sample in the COSMOS field 
and a mid-infrared (MIR) selected sample from the ELAIS 
fields. We will be also revisiting the SDSS/SWIRE quasar 
sample studied in Paper 1 (hereafter quasar sample). 



2.2 X-ray, radio and Mid-IR selected type 2 AGN 

The next sample, hereafter COSMOS sample, consists of 
AGN select ed among X-ray an d radio emitters in the COS- 
MOS field l|Trump et al.ll2007l ). Among the 466 objects for 
which Magellan spectroscopy was carried out, were 135 type 
2 AGN and hybrids (type 2s with a galaxy component; for 

1 http: / /www. mpa-garching.mpg.de/SDSS/DR4/Data/agncatalogue. html 

2 The four IRAC and three MIPS channels' effective wavelengths 
are 3.6, 4.5, 5.8 and 8.0 pm and 24, 70 and 160 um, respectively 
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Table 1. SDSS names, coordinates, redshifts, emission line properties, and stellar masses for the 420 objects overlapping the SWIRE 
fields and the narrow emission line object catalogue from Kauffmann et al (2003). Column called "S" (for Subset) marks the objects that 
have an adequate sampling of their SEDs and will be used in the analysis. The full catalogue is available online. 



SqNr S SDSS name RA Dec z log[OIII] log[OIII/Hb] log[NII/Ha] logM*[M Q ] 



1 




SDSS. 


.J160215. 


.17+535456. 


1 


16 


:02: 


15 


.17 


+53: 


51: 


56. 


1 





.081 


6.2711 


0.341 


-0.0537 


11.065 


2 


X 


SDSS. 


J160258. 


.45+540929. 


1 


16 


:02: 


:58 


.45 


+54: 


09: 


29. 


1 





.065 


5.6622 


0.113 


0.1397 


10.4765 


3 


X 


SDSS. 


J160300. 


.04+541737. 


3 


16 


:03: 


00 


.04 


+54: 


17: 


37. 


3 





.106 


5.862 


0.0946 


-0.0224 


11.1844 


1 


X 


SDSS. 


J160221. 


.35+541105. 





16 


:02: 


21 


.35 


+54: 


11: 


05. 








.066 


6.4168 


-0.0782 


-0.2444 


11.2855 


5 




SDSS. 


.J160006. 


.83+541717. 


2 


16 


:00: 


06 


.83 


+54: 


17: 


17. 


2 





.066 


6.4527 


-0.115 


-0.2069 


10.5041 


6 


X 


SDSS. 


.J160114. 


.70+542034. 


3 


16 


:01: 


14 


.70 


+54: 


20: 


31. 


.3 





.064 


7.7295 


-0.1227 


-0.2971 


10.7284 


7 


X 


SDSS. 


J160035. 


.48+541529. 


2 


16 


:00: 


35 


.48 


+54: 


.15: 


29. 


2 





.064 


6.1645 


0.4281 


-0.1508 


10.7402 


8 




SDSS. 


.J155713. 


.43+543953. 


9 


15 


:57: 


13 


.43 


+54: 


39: 


53. 


9 





.045 


5.9737 


-0.0429 


-0.0453 


10.712 


9 


X 


SDSS. 


.J155813. 


.02+543656. 


1 


15 


:58: 


13 


.02 


+54: 


36: 


56. 


1 





.047 


5.6459 


-0.2685 


-0.2906 


10.4376 


10 




SDSS. 


J155705. 


.98+543902. 





15 


:57; 


05 


.98 


+54: 


39: 


02. 








.045 


6.261 


-0.0665 


0.0228 


10.8279 


11 




SDSS. 


J155708. 


.97+544355. 


.7 


15 


:57: 


08 


.97 


+54: 


,13: 


55. 


.7 





.047 


5.6832 


0.1481 


0.0986 


10.7445 


12 


X 


SDSS. 


J160007. 


.28+545350. 





16 


:00: 


07 


.28 


+54: 


53: 


50. 








.084 


7.264 


-0.3575 


-0.247 


10.5545 


13 


X 


SDSS. 


.J155940. 


.66+545050. 


.8 


15 


:59: 


10 


.66 


+54: 


50: 


50. 


.8 





.093 


6.031 


0.2662 


0.0451 


10.9104 


14 


X 


SDSS. 


.J160218. 


.41+544718. 


5 


16 


:02: 


18 


.11 


+54: 


:47: 


:18. 


.5 





.083 


6.887 


0.0743 


-0.1568 


10.6365 


15 


X 


SDSS. 


J160254. 


.28+545540. 


7 


16 


:02: 


:54 


.28 


+54: 


55: 


10. 


.7 





.120 


6.5346 


0.3394 


0.1132 


11.1944 


16 


X 


SDSS. 


J160124. 


.58+550120. 


1 


16 


:01: 


24 


.58 


+55: 


01: 


20. 


1 





.037 


7.1557 


0.7 


0.0261 


10.6611 


17 


X 


SDSS. 


J160512. 


.63+545550. 


.7 


16 


:05: 


12 


.63 


+54: 


:55: 


50. 


7 





.064 


6.0292 


0.7485 


-0.0783 


9.9791 


IS 


X 


SDSS. 


.J160122. 


.30+552213. 


6 


16 


:01: 


22 


.30 


+55: 


22: 


13. 


6 





.119 


6.3901 


-0.3305 


-0.1667 


10.9273 


19 




SDSS. 


.J161320. 


.00+525516. 


.4 


16 


:13: 


20 


.00 


+52: 


.55: 


16. 


.4 





.166 


8.0556 


0.3297 


-0.1401 


11.2084 


20 


X 


SDSS. 


.J161339. 


.10+531833. 


6 


16 


:13: 


39 


.10 


+53: 


:18: 


33. 


6 





.107 


7.3537 


-0.2625 


-0.2766 


11.1697 



Table 2. SWIRE fluxes of the sample presented in Table [T] The typical errors on the 70 and 160 pm fluxes are 10% of the fluxes. The 
column "AGN" denotes the objects that, after the fitting, were found to have an AGN component and on which all AGN-related analysis 
is based. The full sample is available online. 



SqNr 


S3. 6 My] 


S4.5 


My] 


S5.8 My] 


Ss.o My] 


S24 My] 


S70 My] 


Sieo My] 


AGN 


1 














649.83 ± 25.48 








2 


982.65 ± 4.09 


642.35 


± 


3.17 


463.57 ± 9.7 


786.24 ± 8.26 


522.05 ± 12.92 








3 


1601.7 ± 4.57 


1112.9 


± 


4.39 


760.76 ± 10.71 


602.03 ± 5.61 










4 


4056.2 ± 13.12 


2806.7 


± 


6.52 


2756.1 ± 21.18 


8360.2 ± 14.44 


7091.1 ± 21.25 


115570.0 


574590.0 




5 














4851.0 ± 21.13 


15020.0 






6 


1707.0 ± 6.56 


1213.2 


± 


3.69 


1483.1 ± 17.45 


5734.6 ± 11.03 


9659.0 ± 21.62 


129220.0 


218440.0 




7 




1119.1 


± 


7.12 




514.31 ± 14.09 


407.65 ± 16.67 








8 














1964.0 ± 14.73 








9 




826.02 


± 


4.27 




760.54 ± 12.53 


1099.9 ± 16.91 








10 














5237.6 ± 20.65 


92560.0 


360010.0 




11 














480.31 ± 15.1 








12 


1321.2 ± 2.47 


924.4 


± 


2.53 


925.13 ± 7.79 


4547.6 ± 8.6 


5377.6 ± 20.53 


67840.0 


80250.0 




13 


1202.8 ± 2.89 


825.9 


± 


3.01 


515.12 ± 9.01 


699.14 ± 8.67 


169.76 ± 17.27 








14 


1134.7 ± 2.58 


751.64 


± 


2.74 


559.29 ± 6.58 


1201.9 ± 7.39 


932.88 ± 17.15 








15 


1379.4 ± 3.89 


1008.8 


± 


4.39 


698.5 ± 8.38 


1753.8 ± 8.52 


2398.5 ± 16.42 


1620.0 






16 


2722.7 ± 3.44 


1720.6 


± 


3.18 


1365.3 ± 8.46 


1661.1 ± 7.45 


1232.4 ± 13.53 








17 


512.8 ± 2.4 


318.82 


± 


2.88 


231.09 ± 7.23 


182.6 ± 8.5 


171.23 ± 16.42 








18 


929.66 ± 2.98 


681.77 


± 


3.56 


502.83 ± 8.02 


1214.6 ± 9.39 


1153.0 ± 17.27 


6960.0 


2010.0 


X 


19 














2803.6 ± 19.2 


32840.0 






20 


1653.2 ± 2.66 


1281.6 


± 


3.46 


1343.0 ± 8.09 


8447.7 ± 12.28 


9574.3 ± 22.22 


131310.0 


396400.0 





details on the classification see lTrump et alJl2007T ) with op- 
tical and IRAC counterparts. For the optical/near-IR pho- 
tometry we used the COSMOS fi rst release optical and near- 
IR catalogue lCapak et all (|2007l ). The IRAC and MIPS pho- 
tometry was taken from the S-C OSMOS May 2007 release, 
described in ISanders et alll2007l . From the 135 objects, 92 
have a 24 micron counterpart and they alone will be the 
objects to be studied, as will be discussed in Section [4721 

To construct the SED of these 92 objects we used the 
u*- and i*-band CFHT data, the g + , Vj, r + , and z + Sub- 



aru data, the K s KPNO or CTIO data and the IRAC and 
MIPS S-COSMOS data. Photometry from the MIPS 24um 
deep scans was used whenever available. The Bj band data 
were discarded due to the large overlap with ry + -band that 
would lead to overly redundant information and the CFHT 
i* was used instead of the Subaru i + because some of the 
objects were missing from the Subaru ^ + -band data (even 
though, according to ICapak et al1l2007l the latter observa- 
tions are deeper by two magnitudes) . Tables [3] and [4] show 
the COSMOS names, coordinates, redshifts and type ("q2" 
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0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 
Spectroscopic Redshift 



Figure 1. Redshift histogram of the low luminosity low redshift 
AGN sample selected from Kauffmann et al. 2003. The grey shad- 
ing indicates sources with good photometric coverage allowing 
detailed study. 



for type 2 AGN, "q2e" for hybrid and a question mark next 
to the type d e notes a questionable classification) taken from 
Trump et ail 1 2007 ). Table 2, and the Spitzer fluxes (from 
Sanders et alj 2007 ). respectively, for the first 20 objects of 



the sample. The full tables are available online. 

We also use a list of type 2 AGN taken from the mid-IR 
selected spectroscopic sample of galaxies a nd AGN from the 
15um EL AIS- SWIRE survey, presented in iGruppioni et aD 
(2008). From the 203 spectroscopically identified objects in 
the sample (hereafter EL AIS sample), we selected a total 
of 23 type 2 AGN that had sufficient photometric coverage. 
The ELAIS names and Spitzer (SWIRE) flux es for these 
objects can be found in IGruppioni et al.l ((2008), Table 1. 



2.3 Combined sample properties 

Fig.[5]shows the position of the galaxies on the IRAC colour- 
colour diagram S8.0/S4.5 versus S5.8/S3.6. In the left panel, 
the samples in black, red, blue, and green are the low-z, 
COSMOS, ELAIS and the quasar samples, respectively. The 
points inside open circles denote detections at 70um. The 
different loci formed by the four different samples reflect 
the distinct selection tech niques and hence pr operties of the 
objects. As suggested by ISaiina et al.l (2005) from a study 
based on empirical templates, and as will be demonstrated 
shortly, the objects at the lower left of the plot are dom- 
inated by stellar emission; the objects rising in S8.0/S4.5 
while maintaining a low S5.8/S3.6 (i.e., the bulk of the low- 
z sample with about a third of the COSMOS sample lying 
in the same region) are those dominated by PAH emission 
while the objects with both colours rising, towards the upper 
right corner of the plot, are continuum dominated objects. 
The quasar locus occupies a very well confined region o f 
the colour space, as shown already (e.g. lLacv et al.ll2004 ); 
in their majority and up to a redshift of ~2, they lie on 
a straight line of slope of one, simply indicating the rise 
of the torus emission as a power law. We should note that 
th e position of X-ray select ed AGN has already been shown 
bv lCardamone et al.l l|2008h . that sample however contained 
a much smaller fraction of objects with star formation, and 




0.0 0.5 




5.0 1e42 5.0 

IRAC3.6 Luminosity [erg/sec] 




MIPS24 Luminosity [erg/sec] 



Figure 3. Redshift (upper panel), 3.6 pm luminosity (middle 
panel) and 24 um luminosity (lower panel) histograms for for 
the low-2, COSMOS, ELAIS and quasar samples. Colour coding 
follows Fig. [2] The black arrow indicates the location of the low-2 
sample. 



was plagued by a large fraction of objects with unknown red- 
shifts. Interestingly, their colors tend to populate the lower 
right part of the colour space, which is strikingly empty in 
our larger sample. 

One should keep in mind that the various samples have 
not only different observed properties but also very different 
average redshifts, as seen in the right panel of Fig. [2] A 
histogram of the redshift distributions is shown in the upper 
panel of Fig. [3] The middle and lower panels of Fig. [3] show 
the distributions of the 3.6 and 24 um luminosities of the 
various samples. 
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Table 3. COSMOS names, coordinates, rcdshifts and classification type for the COSMOS type 2 AGN sample, from Trump et al. 2007. 
The full catalogue is available online. 



SqNr 



COSMOS Names 



RA 



Dec 



Type 



1 


COSMOS. 


J100033. 


.08+015414 


.4 


10 


:00 


:33 


.08 


+01 


:54: 


:14 


.4 


1 


.16908 


q2? 


2 


COSMOS. 


J095945. 


.18+023439 


.4 


09 


:59 


: 15 


.18 


+02 


:34 


:39 


.4 





.47717 


q2e 


3 


COSMOS. 


J095943. 


.76+022008 


.0 


09 


:59 


:43 


.76 


+02 


:20 


:08 


.0 





.92998 


q2e 


4 


COSMOS. 


.J100035. 


.25+015726 


.0 


10 


:00 


:35 


.25 


+01 


:57 


:26 


.0 





.84845 


q2e 


5 


COSMOS. 


.J100059. 


.47+022631 


.8 


10 


:00 


:59 


.47 


+02: 


:26 


:31 


.8 





.66944 


q2e 


6 


COSMOS. 


.J100204. 


.84+013356 


.0 


10 


:02 


:04 


.84 


+01 


:33: 


:56 


.0 





.09815 


q2 


7 


COSMOS. 


J100013. 


.73+021221 


.3 


10 


:00 


:13 


.73 


+02 


:12 


:21 


.3 





.18645 


q2 


8 


COSMOS. 


.J095933. 


.79+014906 


.3 


09 


:59 


:33 


.79 


+01 


:49 


:06 


.3 





.13327 


q2 


9 


COSMOS. 


.J100028. 


.69+021745 


.3 


10 


:00 


:28 


.69 


+02 


:17 


:45 


.3 


1 


.03911 


q2? 


10 


COSMOS. 


.J100107. 


.22+020115 


.8 


10 


:01 


:07 


.22 


+02 


:01 


:15 


.8 





.24672 


q2 


11 


COSMOS. 


J100223. 


.62+021132 


.8 


10 


:02 


:23 


.62 


+02: 


:11 


:32 


.8 





.12262 


q2e 


12 


COSMOS. 


.J100211. 


.38+015633 


.2 


10 


:02 


:11 


.37 


+01 


:56 


:33 


.2 





.21752 


q2 


13 


COSMOS. 


.J100243. 


.25+020521 


.9 


10 


:02 


:43 


.25 


+02 


:05 


:21 


.9 





.21398 


q2 


Id 


COSMOS. 


.J100022. 


.94+021312 


.7 


10 


:00 


:22 


.94 


+02 


:13: 


:12 


.7 





.18583 


q2 


15 


COSMOS. 


.J100135 


.22+023109 


.1 


10 


:01 


:35 


.23 


+02 


:31: 


:09 


.1 





.21898 


q2 


16 


COSMOS. 


.J095934. 


.89+015241 


.1 


09 


:59 


:34 


.89 


+01 


:52 


:41 


.1 





.13333 


q2 


17 


COSMOS. 


J100152. 


.90+020113 


.7 


10 


:01 


:52 


.90 


+02 


:01 


:13 


.7 





.22042 


q2 


18 


COSMOS. 


.J100118. 


.97+020822 


.5 


10 


:01 


:18 


.97 


+02 


:08 


:22 


.5 





.16768 


q2 


19 


COSMOS. 


.J100231. 


.87+015926 


.7 


10 


:02 


:31 


.87 


+01 


:59: 


:26 


.7 





.80314 


q2 


20 


COSMOS. 


.J095945. 


.18+015330 


.4 


09 


:59 


:45 


.18 


+01 


:53 


:30 


.4 





.21987 


q2e 



Table 4. COSMOS Spitzer fluxes of the sample presented in Table \3\ The errors on the 70 and 160 pm fluxes are 10% of the fluxes. The 
full sample is available online. 



SqNr 


S3. 6 [/^Jy] 


S4.5 I/^Jy] 


S5.8 


[i"Jy] 


S8.0 \pJy] 


S24 [/^Jy] 


S70 liJ-Jy] 


S160 I/^Jy] 


1 


148.19 ± 0.26 


145.3 ± 0.27 


108.12 


± 


0.74 


680.25 ± 1.55 


2819.9 ± 115.8 


48550.0 




2 


329.74 ± 0.4 


239.16 ± 0.42 


155.24 


± 


0.74 


302.6 ± 1.77 


502.5 ± 11.5 






3 


116.34 ± 0.23 


102.27 ± 0.26 


96.91 


± 


0.77 


125.99 ± 1.47 


562.8 ± 10.8 






1 


108.87 ± 0.22 


114.63 ± 0.27 


75.65 


± 


0.74 


235.19 ± 1.52 


687.8 ± 111.5 






5 


10.37 ± 0.1 


18.13 ± 0.17 


34.93 


± 


0.61 


92.23 ± 1.3 


851.0 ± 108.2 






6 


806.09 ± 0.89 


627.66 ± 0.7 


767.83 


± 


1.38 


5846.26 ± 4.76 


8750.9 ± 170.7 


187800.0 




7 


132.39 ± 0.2 


106.82 ± 0.25 


71.29 


± 


0.62 


262.16 ± 1.44 


3295.0 ± 42.3 






8 


376.01 ± 0.45 


279.86 ± 0.39 


267.35 


± 


0.84 


2175.11 ± 2.45 


3684.0 ± 115.0 


73760.0 


219900.0 


9 


33.02 ± 0.14 


41.25 ± 0.2 


43.11 


± 


0.72 


44.6 ± 1.41 


396.7 ± 112.0 






10 


73.22 ± 0.17 


80.59 ± 0.24 


56.14 


± 


0.67 


300.27 ± 1.58 


517.2 ± 96.7 






11 


365.01 ± 0.4 


285.97 ± 0.38 


255.93 


± 


0.77 


1464.55 ± 1.97 


1957.6 ± 108.3 


55820.0 




12 


80.37 ± 0.17 


75.51 ± 0.21 


48.93 


± 


0.61 


292.38 ± 1.39 


654.9 ± 112.3 






13 


242.32 ± 0.33 


228.27 ± 0.31 


172.33 


± 


0.78 


1058.68 ± 1.69 


2214.2 ± 94.0 






14 


246.44 ± 0.31 


203.64 ± 0.33 


159.9 


± 


0.73 


1273.79 ± 2.03 


6399.4 ± 111.3 


63680.0 




15 


90.39 ± 0.19 


84.12 ± 0.22 


57.78 


± 


0.68 


373.63 ± 1.48 


1189.4 ± 222.1 






16 


196.56 ± 0.27 


164.13 ± 0.3 


136.25 


± 


0.67 


489.73 ± 1.55 


1886.0 ± 108.8 






17 


118.25 ± 0.23 


114.71 ± 0.27 


73.56 


± 


0.71 


454.65 ± 1.62 


1164.5 ± 111.0 






18 


285.06 ± 0.4 


231.32 ± 0.36 


178.53 


± 


0.84 


1378.11 ± 2.13 


2221.6 ± 109.6 


42480.0 




19 


22.08 ± 0.12 


19.33 ± 0.18 


18.82 


± 


0.59 


38.77 ± 1.34 


1843.1 ± 94.6 






20 


81.2 ± 0.17 


72.98 ± 0.2 


42.49 


± 


0.67 


142.2 ± 1.21 


432.0 ± 88.7 







3 OBSERVED AND MODEL SEDS 



3.1 The stellar component 



The observed UV to FIR SED of a galaxy can be decom- 
posed in three distinct components: stars with the bulk of 
their power emitted in the optical and near-IR, hot dust 
mainly heated by UV/optical emission from gas accreting 
onto the central supermassive black hole and whose emission 
peaks somewhere between a few and a few tens of microns, 
and cold dust principally heated by star formation. In the 
present work we consider all three components, which we 
model as follows. 



The stellar component is the sum of Simple Stellar Popu- 
lation (SSP) models of different age, all having a common 
(solar) metallicity. T he set of SSPs i s buil t using the Padova 
evolutionary tracks (|Bertelli et al.l 1 19941 ) , a Salpeter IMF 
with m asses in the range 0.15 - 120 M© and the lJacobv et al.l 
( 1984) library of observed stellar spectra in the optical do- 
main. The extension to the UV and IR range is derived 
from the Kurucz theoretical libraries. Dust emission from 
circumstellar e nvelo pes of AGB stars has been added by 
iBressan et al.l (|l99&f ). The weight in the final spectrum of 
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Figure 2. IRAC colour-colour plot (S8.0/S4.5 versus S5.8/S3.6) illustrating the objects of all the samples detected in the four IRAC 
bands. Left panel: the samples in black, red, and blue are the low-z, COSMOS, and ELAIS, respectively. The quasar sample (green) is 
shown for comparison. The points inside open circles denote detections at 70pm. The right panel shows the same diagram but with the 
colours of the points indicating their redshift. 



each SSPs, as a function of age, is computed according to a 
Schmidt law for the star formation rate: 



SFR(t) 



( Tg-t 



x exp 



Tg-t 
Tg ■ T sf 



(1) 



where Tg is the age of the galaxy (i.e. of the oldest SSP), 
which is assumed to be as old as the age of the universe at 
the galaxy's redshift, and r s / is the duration of the burst 
in units of Tg. Extinction is applied to the final SED by 
assuming a uniform fore ground dust screen w ith a standard 
Galactic extinction law ( Cardelli et al.lfl989T >. The free pa- 
rameters are, therefore, the duration of the initial burst and 
the amount of extinction. 



3.2 The torus component 

T he models for du st emission in AGN are those presented 
m iFritz et al.ll200fj . To summarize, the dust distribution is 
smooth, the torus geometry is that of a "flared disk", the 
dust consists of graphite and silicate grains. The graphite 
grains, with sublimation temperature around 1500 K are re- 
sponsible for the black body like emission in the near-IR 
part of the AGN spectra, creating a minimum in the ob- 
served SEDs between the falling accretion emission and the 
rising torus emission, while the silicate grains are respon- 
sible for the absorption feature at 9.7 micron, seen in the 
IR spectra of all type 2 AGN. The size of the inner torus 
radius, Ri n , depends both on the sublimation temperature 
of the grains , X1500 , and on t he accretion luminosity, L acc , 
according to lBarvainisI (|l987h : 



Rin ~ 1.3 ■ 



(2) 



The dust density within the torus is allowed to vary along 
both the radial and the angular coordinates: 



-7 1 cos (8) I 



(3) 



We give here, as a reminder, the parameter values for 
the discrete grid of torus models that will be used, as pre- 
sented in Paper 1: ft = 0.0, —0.5, —1.0 and 7 = 0.0, 6.0 (see 
Eq-0, optical depth T9.7 = 1.0, 2.0, 3.0, 6.0, 10.0; torus open- 
ing angle O=20°, 40°, 60°; and outer-to-inner radius ratio 
R O ut/Rin=30 and 100. This last set of parameters only holds 
when (i = 0.0, i.e. when the dust density is constant with 
the distance from the centre. For the models with < 0.0, 
Rout is re-calculated to be the distance at which the dust 
density drops at 10% of its value at Ri„. 



3.3 The starburst component 

For the cold dust component, which is the major contribu- 
tor to the bolometric emission at wavelengths longer than 
~ 30 um, we choose six well-studied observational SB tem- 
plates: M82 as a representative of a "typical" SB IR emis- 
sion, Arp 220 as representative of a very extinguished star- 
burst and the templates of NGC1482, NGC4102, NGC5253 
and NGC7714, as intermediate SB templates. 



3.4 Fitting the observed SEDs 

The SED fitting procedure is that of a standard \ 2 minimi- 
sation. As a first step, the UV to near-IR (rest-frame) data 
points are fitted by means of a stellar host component that 
dominates the emission in this range (since the accretion 
emission of the low luminosity AGN of the sample is most 
probably masked under the stellar light). The best fit to the 
observed data points in this wavelength range is found by 
exploring a grid of r s / and E(B — V) values (see Sec. 13 - 1 ft - 
In the next step, the mid-IR points are then fit by a 
torus component. Since all the objects under study are type 
2 AGN and therefore the central source is hidden behind 
the torus, we will focus on the results of the run realised 
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with torus models of high optical depth, T9.7 ^ 1.0 with 
the exception of the low-z sample, for which we will allow 
a full run, including rg.7 < 1.0 torus models, since for these 
low luminosity objects obscuration might also be the result 
of the dust in the host galaxy. In order to properly con- 
strain the part of the SED where the cold dust is domi- 
nant, one would need at least one FIR measurement, in this 
case a 70 and/or 160 pm detection (the mid-Infrared 24 pm 
point alone is probably not enough to distinguish between 
a torus- like and a starburst-like emission). Correspondingly, 
an SB component will only be allowed in the presence of a 
70 (and/or 160) pm datapoint(s), along the lines of Paper 
1. A single deviation from this will be discussed in Section 

EED 

As detailed in Paper 1, the SED fitting leads to the 
computation of a series of physical parameters, namely the 
accretion luminosity, L acc , the IR luminosity, Lir, integrated 
between 1 and 1000 pm, the inner and outer torus radii, Ri n 
(defined in Eq.[5]) and R ou t, the optical depth (or extinction) 
and hydrogen column density along the line of sight, the 
covering factor, CF, i.e. the fraction of the radiating source 
covered by the obscuring material, and the mass of dust, 
Moust- A quantity that was not dealt with in Paper 1 due 
to the nature of the sample (type 1 quasars and hence almost 
complete lack of stellar components) is the stellar mass, M* , 
derived directly by the SSPs, that we will also include in the 
analysis of the low-z sample. 

For the degeneracies related to the various torus models 
parameters and the limitations inhibited in our approach, we 
defer the reader to sections 5.6 and 6 of Paper 1, where all 
these issues were discussed in detail. Dealing with quasars, 
the stellar and SB components were of secondary impor- 
tance. In the present study, however, these two components 
play a more important role. The stellar component (SSPs), 
even though used to fit mainly the UV-to-optical part of 
the SEDs whenever necessary, is not studied any further 
and therefore discussing degeneracies between the relevant 
parameters is well outside the sco pe of this wo r k, but a rel- 
evant discussion can be found in iBerta et al.l (|2004h . The 
only remaining issue is that of the normalisation of the SB 
template, particularly in the absence of a sufficient number 
of data points. After an initial normalisation of the men- 
tioned temlpate to the 70 pm point, an iterative process 
results in the best normalisation factor after the torus emis- 
sion has been excluded. In other words, the SB template is 
finally used to fit the emission that is not accounted for by 
the torus model. This could introduce a bias, maximising 
the estimated AGN and minimising the SB contributions, 
respectively. However, since the iterative process takes into 
account not only the best fit torus model but the best 30 
torus models, we are confident that these deviations will not 
be more than a few percent. 



4 TESTING THE INDIVIDUAL SAMPLES 

4.1 Low-luminosity, low-redshift AGN 

The distribution of the reduced % 2 °f the sample is very 
broad and in general the values are quite large (Fig. [4J plain 
line). The reasons for these high values has been explained in 
detail in Paper I and relate, among other things, to the small 
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Figure 4. Distribution of the minimum x 2 from the analysis of 
the low-2 sample from the "standard" run when no SB component 
was allowed in the absence of 70 (and/or 160) um data (plain 
line) and from the run allowing an additional SB component in 
the presence of a 24 um detection (dashed line). In gray, is the 
final x 2 distribution selecting, for each object, the minimum x 2 
between the two runs. 



values of the photometric errors as well as the crudeness of 
the X-ray-to-optical spectrum of the AGN models. 

Examples of the observed SEDs and the best fit models 
for objects assigned a torus template are shown in Fig. [5] 
The sequence number (SqNr) shown in the plots correspond 
to the SqNr shown in tableQ] The signature of a torus would 
be a large increase of the 24 pm flux with respect to the 8.0 
pm point while the presences of a starburst will usually be 
revealed by a rise in the 8 pm flux due to the PAH features, 
but with a subsequent flattening or slow rise of the SED 
towards the redder (24 pm) wavelength. The absence of any 
component other than the stellar is marked by a monotone 
decrease of the SED after the stellar near-Infrared bump 
(~ 2 — 3 pm ). The best fits for all 388 objects are provided 
as online material. 

In Fig. [6] we also provide three examples of objects for 
which the fitting method completely failed to reproduce the 
observed SEDs, in order to give the reader a flavour of the 
various problems we have encountered. Object with SqNr 
124 (left plot in Fig. [6} has very strong mid-to-far infrared 
emission that none of our SB templates was able to account 
for; object 164 (middle plot) has a stellar-like SED up to 
5.8 pm, but the subsequent rise in flux at 8 pm and drop at 
24 pm can not be reproduced by any combination of torus 
and/or SB templates; and the fitting method simply did not 
work for object Nr 367 (right plot), where both the assigned 
SSPs and the torus model reproduce very poorly the ob- 
served SED. 

Going back to the best fits now, there is a correspon- 
dence between the weighted contribution of the various 
model components and the IRAC colours of the objects. Fig. 
[7] illustrates the position of the objects on the IRAC colour 
diagr am (as in F i g. [2t. In agreement with the predictions 
(e.g. ISaiina et al.l I2005T ) almost all objects assigned a stel- 
lar component alone are gathered at the lower left part of 
the diagram (stars), starbursts (large open circles) form the 
bulk of the sample and AGN (objects with a torus compo- 
nent shown as small filled circles) occupy the space towards 
the upper right part of the diagram. Note however, the large 



8 Hatziminaoglou, Fritz & Jarrett 



10 100 




A [/-inn] 



Figure 5. Best fit SEDs for the first 12 low-z objects with reduced \ 2 < 16.0 that were assigned a torus component. Data points (in the 
rest frame) are shown in red; the stellar model denoted with a green dashed line, the torus with a blue line, and the cold dust component 
with a dashed light green line. 



overlap between SBs and stellar-dominated galaxies, which 
suggests their close physical relationship and fitting redun- 
dancy between these two types of galaxies. 

Drawing a conclusion on the nature of the objects with 
no 70 um detection is actually not straight forward. There 
are cases where the IRAC fluxes are monotonically decreas- 
ing, tracing unmistakably the stellar emission. While in 
other objects the 8 and/or 24 um flux is slowly increas- 
ing (unlike a torus signature), indicating the presence of a 
starburst, yet no FIR is detected for some of these cases. 
Objects with predominantly, starburst as opposed to AGN 
activity and SEDs dominated by emission in the MIR (~ 24 



um) rather than the FIR (A ^ 60 urn) a re actually known 
to exist (see e.g. lEngelbracht et al.ll2008l ). 

The observed SEDs of these objects were poorly repro- 
duced with the combination of a stellar and a torus compo- 
nents alone. As a test, and conversely to the process followed 
in Paper I, we performed an extra run allowing for the use 
of an SB template even in the absence of a 70 pm detec- 
tion. The results of this run should, therefore, influence the 
objects without 70 pm counterparts, representing however 
the majority (2/3) of the sample. The distribution of the 
reduced x 2 f° r this run, shown by a dashed line in Fig. [4] is 
narrower than that of the "standard" run, peaking at \ 2 ~ 
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Figure 6. Examples of failed fits for three objects in the low-z sample. 
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Figure 7. IRAC colour-colour diagram (see Fig. [2]l indicating 
the objects with a stellar component alone (stars), with a torus 
component (small filled circles) and a starburst component (large 
open circles). The symbols in red represent objects with a 70 um 
cold dust component. 



2. This improvement comes mainly from better fitting the 8 
and/or 24 pm points of the objects that were only assigned a 
stellar component, by an additional SB component. In order 
to combine the results of the two runs, we will simply select 
the best \ 2 from the two runs, shown in Fig. [4] with a gray 
histogram. For consistency with the study in Paper 1 we will 
consider "good" all fits with \ 2 < 16, whose limit includes 
~70% of the sample. Note, however, that we run the same 
analysis on the entire sample (388 objects) and found none 
of the results differ in any way if all fits are included. 

The sample of objects harbouring an AGN according to 
the present analysis, shows very specific trends: the major- 
ity of objects (87%) were assigned very high optical depth 
tori (T9.7 = 10.0), and 90% were assigned tori models with 
decreasing dust density from the centre (/? < 0.0 in Eq. [3}- 
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Figure 8. Histogram of the fraction of IR luminosity attributed 
to the AGN, for objects assigned a torus component (open his- 
togram) and those with an extra SB component (gray region). 



In fact, there is a clear preference for models with (3 — —0.5, 
with 72% of the objects assigned such a torus component. 
Almost 85% of the sample have small outer to inner radius 
ratio (Rout/Rin = 30). The combined result is very small 
tori with inner radii of less than 1.4 pc ( 0.185 pc average) 
and outer radii reaching in a few cases 100 pc. These average 
values are ~2 times smaller than the average values for Ri„ 
derived for the quasars and reflect the lower average accre- 
tion luminosity, L aC c, of the low-z sample. In fact, the lower 
(by a factor of ~25) (L a cc) of this sample with respect to 
the quasars is not only due to the different redshift ranges 
covered by the two samples but also to the very nature of 
the objects. 

The contribution of the AGN to the IR luminosity, for 
the 75 "good" fits that have been assigned a torus compo- 
nent, is shown in Fig. |8j with the grayed-out region corre- 
sponding to the 20 objects with an extra SB component. 
Due to the nature of the objects the AGN fraction tends 
to be very low compared to the values derived from the 
SWIRE/SDSS quasar sample of Paper 1. 

All 91 objects with a 70 pm detection, for which the fits 
are optimally constrained, were assigned a starburst compo- 
nent; of that total, 20 of them show evidence for the presence 



10 Hatziminaoglou, Fritz & Jarrett 



0.6 

0.4 

0.2 

-0.0 

■0.2 

-0.4 

-0.6 
-0.8 

















































® 


> 




AG 

) 


N 
















■> 














— « 




® 


















- TT 








< 


^ — S 








C 


® 


o 




) 






■ 

















O 


LINERS 














— » 










H 


















Star 


burst 


°c 





































8.5 

8 

7.5 

7 o 
m 

O 

6.5 = 



5.5 



-0.5 -0.4 -0.3 -0.2 -0.1 -0.0 0.1 
log(NII/Ha) 



Figure 9. Emission line ratios as a function of extinction cor- 
rected [Om] line luminosity, for those objects with a 70 um coun- 
terpart. Large open and small filled circles denote SB and torus 
component, respectively. The dashed line shows the usual - and 
more conservative - AGN selection criterion suggested by Kew- 
ley et al. (2001), according to which the AGN should occupy the 
upper right part of the graph. 



Table 5. The average values for the [Om] luminosity and the 
log([Oin]/H/3) ratio for the objects with an AGN component, an 
additional SB component and an SB component alone. 



Component 


# 


log([Oin]) 


log([Om]/H/3) 


AGN 


20 


8.18 ± 0.60 


0.299 ± 0.42 


AGN + SB 


91 


7.38 ± 0.75 


0.006 ± 0.37 


SB alone 


71 


7.16 ± 0.64 


-0.077 ± 0.31 



of a torus. Fig. [9] depicts the relative intensities of emission 
lines (log([Om]/H/3) versus log([Nll]/Ha), as a function of 
the extinction corrected [Om] luminosity. Six of nine objects 
that occupy the space commonly agreed to be populated by 
Seyfert 2 galaxies (upper right corner in Fig. |5J) have indeed 
been assigned an AGN component. Another four lie very 
close to the dividing lines, and a total of six in the LIN- 
ERs' region (lower right part). The majority of the objects 
were assigned an SB component alone and w ould have not 
been p art of the AGN catalogue at all, had theE ewlev et al.l 
|200ll ) AGN selection criterion been applied (shown here in 
a dashed line). 

The average values for the [Om] luminosity and the 
log([Oni]/H/3) ratio for the objects with an AGN compo- 
nent, an additional SB component and an SB component 
alone are given in table [S] Both quantities have larger val- 
ues when an AGN component is present and lower when 
an SB component alone is responsible for the cold(er) dust 
emission. 

Our analysis is only able to identify a torus compo- 
nent in about one-third of the total low-z sample (75 [142 
if no cut in the \ 2 1S imposed] out of 266 [388]). Even 
though evidence of ch anging behaviour of (clumpy) tori 
l|Honig fc Beckertj|2007t) or even complete disappearance of 
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Figure 10. Covering Factor distribution for the low-z sample. 



the torus |Elitzur fc Shlosman 2006) towards the lower lumi- 
nosity regimes (L < 10 42 erg/sec) has been reported, these 
effects are not what we observe here, as the lowest accretion 
luminosity estimated for objects with a torus component is 
of the order of 10 44 erg/sec, seen in the first plot of Fig. 1181 
We therefore conclude that we either are unable to detect 
the signature o f AGN with L acc < 10 4 4 erg /sec or that the 
criteria used bv lKauffmann et all (|2003D to select AGN were 
too generous and many objects with no nuclear activity may 
have made it into the suggested AGN sample. 

Another issue to consider is the distribution of the re- 
sulting covering factor, shown in Fig. 1101 Were the com- 
puted values for CF correct, they would translate to a ratio 
of type2:typel objects of about 1:1 in the local and low-z 
universe, a much lower value than that observed. 

This discrepancy could again be seen as a result of the 
AGN candidate selection and the somewhat sparse samples 
of the observed SEDs, the combination of which does not 
always allow the fitting code to give robust results. 

Allowing for torus models with T9.7 < 1.0 does not have 
any effect on the results, as only 18 out of the 388 objects 
where better fit by a model including such a torus compo- 
nent, with a single one among them assigned a \ 2 < 16.0. 



4-1.1 Properties of the host galaxies 

The SSPs used to reproduce the blue part of the SEDs pro- 
vide information about the stellar mass of the host galaxy. 
To compute the stellar mass we exploit the fact that SSP 
spectra are given in luminosity units per one solar mass. 
We also account for the fact that a certain percentage of 
stars, depending on the SSP's age, has evolved and is not 
shining anymore, so the mass value we provide is actually 
"luminous mass", Fig. [TT] shows a comparison of the stellar 
mass computed in t his work and that provid ed in the narrow 
line AGN catalogue [Kauffma nn et alj l|2003l ). The line shows 
the linear correlation between the two quantities, M^ ere — 
1.59 + (1.06Mk). Note that the systematic differences in 
the mass values are the result of the different IMF used by 
iKauffmann et alj (|2003l ), a Kroupa (2001) IMF, which dif- 
fers from ours both in the slope and in the mass limits. 

We now compute t he star formation rate (SFR) accord- 
ing to lKennicutl l| 19981 ) as 
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Figure 11. Comparison of the stellar mass derived from this 
work based on SED fitting and that presented in Kauffmann et 
al. 2003. 
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Figure 12. Star formation rates for all objects with an SB com- 
ponent and x 2 < 16.0 (black histogram) and those without an 
AGN component (gray histogram). 



SFR[M Q /yr] = 4.5 x 1(T 44 x Lfir (4) 

where Lfir is the infrared luminosity of the SB compo- 
nent, integrated between 8 and 1000 micron. Fig. 1121 shows 
the distribution of the SFRs for all objects with an SB com- 
ponent and x 2 < 16-0 (black histogram), with the grayed 
region corresponding to the objects without an AGN com- 
ponent. The average SFR of the entire sample without im- 
posing a cut on the \ 2 values is of 2.8 Mq/yv, and drops to 
2.2 when only objects with x 2 < 16-0 are consider, and to 
1.85 when objects with an AGN component are excluded. 

Note that all objects with a torus component are gath- 
ered in the part of the histogram with the higher SRFs. We 
will get back to the discussion on SFRs in Section \5. 41 
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Figure 13. Histogram of the minimum \ 2 from the analysis of 
the COSMOS (solid line) and the ELAIS (dashed line) samples. 



4.2 Type 2 AGN and hybrids from the COSMOS 
sample and Mid-IR selected ELAIS AGN 

The spectroscopic sample of X-ray emitte r s sele cted in the 
COSMOS field presented in iTrump et all l|2007l ) comprises 
80 type 2 AGN, 47 type 2 AGN and red galaxy hybrid and 8 
type 2 AGN whose classification is not 100% secure. Ninety 
two out of these 135 objects have a 24 micron counterpart 
and, as already mentioned in Section [2] we will be concen- 
trating on this sub-sample alone. This was not a constraint 
applied to the low redshift sample, however, due to the aver- 
age (low) redshift of this sample as well as its nature (see Fig. 
[3]). The SEDs are dominated by the stellar light all the way 
to the IRAC3 band (5.8 um) and even though the presence of 
a torus is somewhat revealed from the lower Ss.o/Ss.g colour 
(with respect to objects without a torus component), a sin- 
gle data point at 8.0 um (IRAC4) is not enough to constrain 
the torus properties. For the record, 36 of the remaining ob- 
jects (those without a 24 um detection) were indeed assigned 
a torus component while six others were better represented 
by SSPs alone. 

Fig. [13] shows the minimum \ 2 distribution for the sub- 
sample of 92 objects whose SEDs extend at least up to 24 
um, (solid line), excluding two objects with minimum \ 2 > 
90 that are obviously erroneous and will not be taken into 
account in what follows. Again for consistency, we will adopt 
a cut-off x 2 °f 16, keeping in mind though that including the 
rest of the fits will not influence the results in any significant 
way. 

Examples of SED fits of objects assigned an AGN com- 
ponent are shown in Figs. [14] and [15] for the COSMOS 
and ELAIS samples, respectively. The fits for all objects are 
available as online material. 

Note that in the upper and lower left plots in 
Fig. [H (objects COSMOS_J100137.69+002844.1 and COS- 
MOS_J100249.59+013520.2, respectively), the observed 24 
um points lie well above the model (in both cases consist- 
ing only of stars and a torus). Such cases indicate the pos- 
sible existence of an additional SB component, which we 
did not use for these fits because of the lack of observed 
70 and /or 160 um data points, as discussed in detail in 
Section 13.41 Similarly, for the lower right object in Fig. 1151 
(ELAIS_J003718-421924), the fit fails to reproduce the 160 
um point, despite the simultaneous existence of a torus and 
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Figure 14. The best fits for the first 12 COSMOS objects with reduced x 2 < 16.0 that were assigned a torus component. Data points 
(in the rest frame) are shown in red; the stellar model denoted with a green dashed line, the torus with a blue line, and the cold dust 
component with a dashed light green line. 



a SB components. This case shows that either the SB tem- 
plates used are not always adequate or, most probably, that 
we are lacking a diffuse component that would account for 
the cold dust in the host galaxy. 

The fitting results show a clear preference for models 
with dust density decreasing with the distance from the cen- 
tre, with 95% of the objects finding a better fit with j3 < 0.0 
models (75% with (3 = —0/5). A further dependence of the 
dust spatial distribution on the altitude from the equator 
was found, with 87% of the objects having a best fit model 
with 7 = 6.0. The resulting tori sizes do not depend on the 
type of object (AGN or hybrid) and the estimated average 



inner tori radii, Ri„, are of 0.85±0.50 and 0.91±0.58 pc, 
for AGN and hybrids, respectively, and with about 60% of 
the objects better matching torus models with small sizes 
(Rout / Rin =30) . More than half of the objects were found to 
have a large covering factor (> 80%). A comparative study 
on the covering factor derived from the various samples will 
follow shortly. 

Once again, the position of the objects in the IRAC 
colour diagram (Fig. [2j and the contribution of the various 
components in their global SEDs are in very good agree- 
ment. Fig. 1 161 shows the IRAC colour diagram for the COS- 
MOS (circles) and ELAIS (squares) samples. 
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Figure 15. The best fits for the first 12 ELAIS AGN with reduced x 2 < 16.0 that were assigned a torus component. As in Figs. 5 and 
13, data points are shown in red; the stellar model is marked with a green dashed line, the torus with a blue line, and the cold dust 
component with a dashed light green line. 



The majority of objects with a 70 micron detection and 
an assigned SB template gather in the PAH-dominated lo- 
cus (upper left part of the diagram) while the AGN con- 
tribution to the IR luminosity increases as objects move 
along the continuum-dominated locus, becoming redder in 
the S5.8/1S3.6 axis. 

From the 23 objects composing the MIR selected type 
2 ELAIS AGN sample, 21 were found to have a torus com- 
ponent, all of which with T9.7 ^ 6.0, 7=6.0 and f3 < 0.0, 
among which 15 found a better fit with /3 — —0.5 templates. 
The average estimated sizes of the tori were very similar to 



those of the COSMOS sample. The reduced x 2 distribution 
for the fits is shown in Fig. \T3[m dashed lines. 

The AGN contribution to the IR luminosity for both 
samples is shown in Fig. 1171 with the grayed region repre- 
senting the objects with an additional SB component. The 
plain (dashed) line and grayed (black) region correspond to 
the COSMOS (ELAIS) sample. Again, the existence of such 
a component can be, in general, well constrained only in the 
presence of data points longward A = 24um and this is valid 
for only 25% and 50% of COSMOS and ELAIS samples, 
respectively. 
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Figure 16. IRAC colour diagram for the COSMOS (circles) and 
ELAIS (squares) AGN. Small filled symbols denote an AGN com- 
ponent, large open symbols an SB component, for objects with a 
good x 2 - Small open symbols mark the objects whose fits were 
not acceptable. 
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Figure 17. Histogram of the fraction of IR, luminosity attributed 
to the AGN, for objects from the COSMOS and ELAIS samples.. 



5 COMBINED AGN PROPERTIES 

We will now attempt to put all the results from the various 
AGN samples together and construct a common description 
of their dust properties. To this end, we will first present a 
summary of the results of Paper 1, that will also be used 
in this combined study. The low- 2 sample will be excluded 
from this analysis, because of the unconfirmed (stellar vs. 
starburst) nature of the objects. They will be taken into 
account, however, when discussing the star formation, as 
the results on this issue seem to be more robust, and only 
marginally affected by the presence of tori, as suggested from 
Fig. El 



5.1 Type 1 SWIRE/SDSS quasars 

Paper 1 focused on the study of 278 SDSS/SWIRE quasars 
and the properties of dust surrounding them. Comparing 
the results between the "restricted" run, where only high 
T9.7 models were allowed and the "full" run, where optical 
depths as low as 0.1 were allowed, the hypothesis of the ex- 
istence of low optical depth tori could not be ruled out. In 
fact, the majority of objects found a better match of their 
SEDs with low T9.7 models (rg.r < 1). The computed aver- 
age inner radius was of ~2.6 pc with a tendency for higher 
Rout/Rini (=100) for the full run (59%) with respect to the 
restricted run (46%), as reported in Paper 1. The covering 
factor depended of course on the choice of the optical depth 
but had a relatively fiat distribution in both cases, taking 
values from as low as < 0.1 to as high as 0.95. As for the 
dust density, there was no clear preference between 7 = 0.0 
and 7 — 6.0 models, in any of the runs, while /3 < 0.0 mod- 
els were always favoured. The accretion and IR luminosities 
were the two better constrained quantities and were almost 
independent on the choice of T9.7. The subsample of 70 ob- 
jects with 70 um detections enabled an estimation of the 
contribution of a starburst to the total IR luminosity, indi- 
cating contributions of up to 80%. 



5.2 Comparative study of type 1 & 2 AGN 

Section 12.31 described the observed properties of the com- 
bined AGN sample. Although the sample is not complete, 
the combined study of the various AGN sub-samples com- 
posing it, will allow for comparative conclusions. 

Fig. 1 181 shows the distribution of the various parameters 
for the different samples. From top to bottom and left to 
right, and following the colour-coding introduced in the left 
panel of Fig. [2] we present: the accretion luminosity; and 
the IR luminosity attributed to the torus component (i.e. 
excluding the stellar emission coming either from the SSPs 
and/or the SB components); the covering factor; the viewing 
angle; the hydrogen column density along the line of site; 
and the mass of dust confined inside the torus. The open 
histogram corresponds to the "restricted" run, where only 
high optical depth tori models (T9.7 ^ 1.0) were allowed 
(Paper 1). 

The combined quasar and type 2 AGN samples cover 
four orders of magnitude with the quasar L a cc dis- 

tribution peaking at luminosities of about an order of mag- 
nitude brighter than the type 2s (see Fig. 1181 - upper left). 
This observational bias affects, for example, the fraction of 
the IR luminosity attributed to the torus component (Fig. 
[T8l - upper right). 

The viewing angle as derived from the fits, shown in 
the right column, second row, is consistent with the Unified 
Scheme postulating that the differences between type 1 and 
type 2 objects are a line-of-sight effect, with the type 2 ob- 
jects seen through the obscuring material. Almost all type 
2 objects are seen in small viewing angles, i.e. close to the 
equator and through large amounts of dust, while type 1 ob- 
jects are seen in large viewing angles especially when only 
tori with high optical depths are considered (open green 
histogram). Note that low optical depth tori would result 
in low obscur a tion t ype 2 AGN when seen edge-on. Fig. 4 
in iFritz et al.l (|2006l ). e.g., shows the emission of an AGN 
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Figure 18. From top to bottom and left to right: accretion luminosity, L aC c, and infrared luminosity attributed to the AGN component, 
covering factor, viewing angle, hydrogen column density, and mass of dust, for the various samples. Green filled and empty histograms 
show the distributions for the quasar "full" and "restricted" runs, respectively; red open full-line and blue dashed-linc histograms show 
the results for the COSMOS and ELAIS samples, respectively. The low-z sample has been excluded. 



through a Tg.7 = 0.1 torus for all lines of sights. While the 
object is clearly a type 1 when seen face-on, the UV-to- 
optical nuclear emission can be absorbed by several orders 
of magnitude when seen through the dust, depending on the 
wavelength and the line of sight. 

The majority of type 2 objects are seen through high 
column densities (Nh 10 23 cm -2 ), while the few quasars 
that are seen through the torus, are actually seen through 
low optical depth and lower than the type 2 column den- 
sity gas (Fig. 1181 lower left panel). Also consistent with the 
Unified Scheme is the distribution of the mass of dust (right 
column, fourth row in Fig. I18[) . If the mass of dust can be 
considered a linear tracer of gas, this implies that the gas 



reservoirs in the vicinity of type 1 and type 2 objects are 
comparable. 

Note that the properties of the radio and X-ray (COS- 
MOS) and MIR-selected (ELAIS) AGN are remarkably sim- 
ilar with the only notable difference being the slightly higher 
IR luminosities of the torus components of the ELAIS sam- 
ple (upper right plot in Fig. I18p . simply reflecting the selec- 
tion effects. 



5.3 The covering factor 

The distribution of the values for the covering factor (seen 
in the left column, second row in Fig. I18|l indicates that tori 
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with both high and low covering factors are possible in both 
types of objects. The COSMOS and ELAIS samples show a 
stronger tendency for higher values. 

The estimated type2:typel ratios based on the CF dis- 
tribution of the type2 objects is ~ 2-2. 5:1; the same ratio is 
found when considering the restricted run of the quasar sam- 
ple. The full type 1 run, however, suggests a ratio of 1:1, with 
a mean value of the covering factor of ~0.47, very c l ose to 
the value of 0.4 computed by iRowan-Robinson et ail (|2008T ) 
for the entire SWIRE quasar population. This low value, 
seemingly in conflict with the Unified Scheme, actually re- 
flects the difference in the accretion luminosity between the 
type 1 and type 2 samples (Fig. 1180 . already discussed in 
Section 1531 

From the type 1 quasars with redshifts below ~ 0.8 we 
can measure the flux at 5100 A, which we can then scale with 
the accretion luminosity, L acc , as seen in the left panel of 
Fig. 1191 Assuming this relation holds for all redshifts as well 
as for type 2 objects, we can translate L acc to Lsioo- We can 
then compute the fraction of obscured AGN in bins of Z/5100, 
shown in the right panel of Fig. \T§\ The black circles and 
open squares make use of the results of the full and restricted 
runs, respectively. The full and two dashed lines show the 
results an d uncertainties ( due to bolometric corrections) ob- 
tained bv lMaiolino et all (|2007l ) who conducted a study of 
25 high luminosity high redshift (2 < z < 3.5) quasars with 
IRS spectroscopy. The observed trend, i.e. the decreasing 
average torus covering factor with increasing AGN luminos- 
ity, is_2n_sjrp_pOTt_ofthe receding torus paradigm, suggested 
bv lLawrence A.lll99ll . 

Despite the crude conversion of L aC c to L5100 our re- 
sults making use of the full run are in very good agree- 
ment (i.e. within the uncertainties) with those derived by 
iMaiolino et al.l (120071 ) based on totally independent method 
and sample, while the results based on the restricted run, 
even though still consistent (excluding the brightest two 
bins) show a quite different trend. We can therefore still 
not conclude on the relative occurence of low T9.7 tori with 
respect to high Tg.7 tori but, based on the above, we can still 
provide some evidence for their existence . 

Previous studies (e.g. IMaiolino et al.ll2007l ) assume that 
the ratio between the primary AGN radiation (here mea- 
sured by Lacc) and the thermal infrared emission attributed 
to the AGN is a direct indicator of the torus covering factor. 
Our findings confirm this assumption, as shown in Fig. [20] 
where the accretion-to-IR AGN luminosity of each individ- 
ual object of each sample is plotted against the computed 
covering factor. The slope of the correlation is very similar 
for type 1 and type 2 samples. Type 2 objects, however, 
have considerably larger L aC c/ Lir(AGN) ratios than the 
quasars. This "jump" in the values is due to the different 
viewing angles and occurs right when the line of sight inter- 
cepts the obscuring torus. In order to illustrate this effect, 
we show the L acc / 'Lir(AGN) ratios as a function of the 
viewing angle for the extreme models of both type 1 and 2 
objects (higher and lower red points; higher and lower green 
points in Fig. [20) in Fig. [21] 

5.4 Starburst activity in AGN 

The study of objects with at least one data point at wave- 
lengths longward of A = 24 um shows that the IR emission 
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Figure 20. Accretion-to-IR, luminosity attributed to the AGN 
component as a function of the covering factor. The green solid 
and dashed lines obtained as linear fits corresponding to the full 
and restricted runs, respectively, for type 1 objects, while blue 
and red lines correspond to ELAIS and COSMOS type 2 AGN, 
respectively. 

of AGN can not be attributed to the torus component alone. 
In the vast majority of cases, and in order to reproduce the 
70 (and, whenever available, 160) um points, an additional 
SB component was necessary. Furthermore, in all those cases 
the contribution of the AGN emission in the infrared is typ- 
ically smaller than 50%, as seen in Figs. and [17] for the 
type 2 objects, and a bit higher in type 1 quasars (Paper 1, 
Fig. 12). 

Fig. [22] shows the estimated SFRs for the various sam- 
ples (with the same colour coding as in Fig. [2]) computed 
from Eq. 3] for all objects with an SB and an AGN com- 
ponents. The average SFR for the low redshift sample is 
about 10 M Q /yr, while those computed for the COSMOS 
and ELAIS samples are of ~40 Mq/jt but with a larger 
span. That of type 1 quasars (shown in green), however, has 
an average SFR of 115 M©/yr, with some cases reaching 
SFRs of the order of 10 3 Mq/ jt. These values are well in 
agreem ent with the findings of ISerieant fc Hatziminaogloul 
(2009), a study based on stacking analysis of a variety of 
quasar samples. These values do not necessarily imply that 
star formation activity is stronger in type 1 quasars in gen- 
eral. More likely, it reflects the observational biases of the 
samples, both in terms of average redshifts and luminosities 
(Figs. [3]and [18]). 

The starburst-to-AGN luminosity ratio for the com- 
bined sample shows a slight tendency to decrease with in- 
creasing accretion luminosity, as seen in Fig. 1231 This is in 
apparent agree ment with previous wo rks on type 1 quasars 
alone (again by IMaiolino et alj l2007h , and if confirmed for 
the entire sample, would imply that type 1 and type 2 ob- 
jects behave in the same way in this respect. On the other 
hand, we can not assert beyond doubt that there there is no 
implicit dependency on the redshift. In fact, examining the 
quasar sample alone (open symbols) we would be tempted 



Dust ton in AGN 17 




Figure 19. Left panel: L5100 measured on the quasar SDSS spectra for objects with redshift lower than ~0.8, versus L aC c, computed 
from the SED fitting. The linear correlation, shown in a solid line, is given by Zog(Lsioo) = 0.692 X log(L acc ) — 1.5767. Right panel: 
fraction of obscured AGN in bins of Lsioo- Filled circles (open squares) represent the results of the full (restricted) run. The solid 
line represents the fraction of obscured AGN as found by Maiolino et al. (2007), while the dashed lines trace the uncertainties due to 
bolometric corrections. 
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Figure 21. Accretion-to-IR luminosity attributed to the torus 
as a function of viewing angle for the type 1 (green dashed 
lines) and 2 (red solid lines) objects with the highest and lowest 
Lace/ Lir(AGN) in Fig. 19. The differences in this fraction are 
mainly attributable to dust self absorption which is more efficient 
for type 2 line of sights. 



to claim an increase with L acc , but which is likely attributed 
to the L - z degeneracy. 




5.0 1e01 5.0 1e02 5.0 1e03 5.0 1e04 5 



SFR [Msun/year] 

Figure 22. SFRs estimated from the Kennicutt (1998) formalism 
for the various samples. 

6 DISCUSSION 

This work focuses on the IR properties type-2 AGN, the 
type 1 issues having already been discussed in a previous 
paper (|Hatziminaoglou et al.l [20081 ). The collection of data 
used for this paper comes from: 1) a low-redshift (0.02 < z < 
0.2) AGN sample, selected via standard emission-lines ratios 
criteria; 2) a higher redshift (0.1 < z < 1.4) sample from the 
COSMOS and ELAIS fields, selected with X-ray and Mid- 
Infrared criteria. The observed SEDs were reproduced by 
means of a three component model, including emission from 
stars, hot dust from AGN and warm dust from starburst 
activity, and acceptable fits were obtained for the majority 
of the objects. 

Despite the fact that the samples we analyzed were not 
chosen homogeneously, a comparison of the behaviour be- 
tween the physical properties of the type 1 and type 2 AGN, 
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Figure 23. Starburst-to-AGN infrared luminosity ratio as a func- 
tion of the accretion luminosity, L aC c, and the redshift. Filled and 
open symbols depict the type 2 and type 1 objects, respectively. 



shows an overall agreement. In Paper I we found the cov- 
ering factor to have a broad trend towards lower values as 
the accretion luminosity, i.e. the power of the central engine, 
increased with small, on-average, covering factors. Applying 
the Lace versus L5100 correlation derived for quasars in Pa- 
per I, we derived the flux at 5100 A for the type 2 samples 
and checked the correlation between the optical luminosity 
and the covering factor, which in turn, can be converted 
into a typel-to-type2 fraction. The results provided by the 
SED fitting analysis, show a very good agreement with the 
relation presented in M aiolino et al. (2007) , who find a cor- 
relation between the fraction of obscured AGN as a function 
of the optical luminosity (i.e. at 5100 A). Another indication 
of the dusty torus was thought to be ratio between the ac- 
cretion luminosity and the torus infrared luminosity. In fact, 
the latter is just reprocessed radiation which scales linearly 
with the primary source: what makes the difference is the 
fraction of "heating radiation" which is intercepted by the 
dust, which again depends on the covering factor. We find 
a very similar trend for these quantities, in both type 1 and 
type 2 objects, also showing how differences in the amount of 
obscuration are very well explained by dust self-absorption, 
i.e. thermal dust emission absorbed by dust iteself. 

For the low-z sample, we found no evidence of emis- 
sion from an AGN component for ~ 70% of the objects. 
Although we cannot rule out the absence of an AGN from 
these sources, we can set an upper limit to its luminosity 
since it is not observed at mid-infrared wavelengths where 
its emission is the strongest as compared to both that of 
the stellar and of the starburst-heated dust. This sample 
suffers from high contamination from the stellar emission 
component, even at mid-infrared wavelengths where (on the 
contrary) it's the AGN component that usually dominates 
when present. The question to ask, therefore, is what can 
be really constrained with three or four data points. One 
of the most reliable quantities should be the optical depth, 
since low values would make the torus emission stand over 
the stellar (SB or stars) in the MIR. In this respect, the 



1ow-jz sample may be too generous in selecting really active 
AGN. There is, sometimes, no evidence for a MIR excess at 
all with respect, for example, not only to a "normal" star- 
burst emission, but also to a passively evolving galaxy. The 
absence of any evidence for a hot dust emission in the MIR, 
which is the place where the stellar component is less impor- 
tant and the AGN contribution is increasingly brighter, will 
remove the AGN component from the fits even though this 
could also be turned into an upper limit. In the cases where 
there is a MIR excess with respect to a pure stellar emis- 
sion, we explore two possibilities: torus and PAH emission. 
In some cases a better fit was obtained adding PAH -i.e. 
starburst component alone- to the model inste ad of AGN 
which should instead dilute the PAH emission (|Lutz et al.l 
Il998h . 

Comparison between clumpy an d smooth tori models 
(e.g. iDullemond fc van Bemmelll2005r ) indicate that globally 
the SEDs produced by the two models are quite similar, but 
with some details characteristic for one or the other model: 
the silicate feature observed in absorption in objects seen 
edge-on is shallower for clumpy models, the average near- 
IR flux is weaker in smooth models and the clumpy models 
tend to produce slightly wider SEDs at certain inclinations. 
Furthermore, clumpy models c an produce very sm all tori 
sizes with R ou t/Rin ~ 5 - 10 (|Nenkova et alj l2008h . while 
still producing a broad MIR emission. Subsequently, the se- 
lection of a smooth torus for the present study might have 
resulted in overestimated tori sizes without, however, jeop- 
ardising the estimates on the properties related to the IR 
emission or our conclusions on the Unification Scheme. SED 
fitting may not be a sensitive enough method to distinguish 
between the differences introduced by the two approaches, 
because of the width of the filters and the scarce sampling of 
the observed SEDs, but also because the various characteris- 
tics of the torus component can be altered or diluted by, for 
example, the presence of a starburst component. Notwith- 
standing these limitations, SED fitting is still the best tool 
available for extracting the maximum information from large 
photometric AGN samples and is now proving to be a power- 
ful technique in relating the dust properties to the accretion 
properties as well as the properties of the larger host galaxy. 

Because AGN are of order hundred times less nu- 
merous than galaxies and also prone to selection biases, 
such AGN studies are only possible with the advent of 
multi-wavelength surveys that probe large volumes and al- 
low the construction of well-samples SEDs. Even though 
we cannot address all issues related to dusty tori, our 
understanding of their properties has been greatly im- 
proved over the last decade thanks to the Spitzer Ob- 
servatory. This space infrared telescope allowed, among 
other things, the construction of SEDs of hundreds of 
AGN of all t ypes thanks to the IRAC and MIPS pho- 



tic 

tome t rv (e.g.lPranceschini et a l. 2005; Hatzi minaoglou et al.l 



120081 : iRichards et al.l l200rj ; IPolletta et al.l |2008f >. Spitzer 
also opened the doors to detailed and coherent studies of 
the interplay between AGN and starburst activity (e.g. 
iHernan-Caballero et al.ll2009l and references therein) , paving 
the way for Herschel. 

In fact, Herschel with it s two cameras /medium r eso- 
lution s pectrographs (PAC S; iPoglitsch fc Altieril [20091 and 
SPIRE: [Griffin et al.l [2009h and a very high resolut ion het- 
erodyne spectrometer fHIFI: |pe Graauw et alj|2009l ) will be 
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the first space facility to completely cover the range between 
60 and 670, where the the bulk of energy is emitted in the 
Universe, allowing for a more detailed sampling of the ob- 
served FIR SEDs and the cold dust emission. Large parts of 
the Herschel key science will focus on the formation and evo- 
lution of galaxies and the studies of star formation. HerMES, 
the Herschel Multi-tiered extragalactic Survey, is the largest 
project that will be conducted by Herschel, with dedicated 
~900 hours that will map over 70 square degrees including 
most of the SWIRE fields (where most of the IR photom- 
etry of all the objects in this study comes from) as well 
as other extragalactic survey fields covered by several mis- 
sions in all wavelengths (e.g. COSMOS or the Groth Strip), 
and will address among others the i s sue of AGN and star- 
burst connection ijGriffin et al.l 120061 : lHatziminaoglou et al.l 
2007). Additionally with Herschel, the Astrophysical Ter- 
ahertz Large Area Survey (ATLAS) will cover about 500 
square degrees and will observe all SDSS quasars with red- 
shift z < 0.2 as well as all the brightest FIR SDSS quasars 
(with luminosities 10 times larger than the mean), amount- 
ing to ~330 individual detections of qua sars in the area 
covered l|Serieant fe Hatziminaogloul I2009T ) (a factor of ~5 
more than all the SDSS quasars with 70 pm detections). 
These studies will allow us to significantly improve our un- 
derstanding of the AGN phenomenon, disentangle the con- 
tribution of cold and hot dust emission in their SEDs and 
study the concomitant occurence of nuclear activity and star 
formation. 
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